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Abstract. We study generation of magnetic fields, involving large spatial scales, 
by convective plan-forms in a horizontal layer. Magnetic modes and their growth 
rates are expanded in power series in the scale ratio, and the magnetic eddy diffusiv- 
ity (MED) tensor is derived for flows, symmetric about the vertical axis in a layer. 
For convective rolls we demonstrate that MED is never below molecular magnetic 
diffusivity. For cell patterns possessing the symmetries of a rectangle, critical values 
of molecular magnetic diffusivity for the onset of small- and large-scale magnetic field 
generation are the same. No instances of negative MED in hexagonal cells have been 
detected. A family of plan-forms has been found numerically, where MED is negative 
for molecular magnetic diffusivity over the threshold for the onset of small-scale mag- 
netic field generation. However, the region in the parameter space, where large-scale 
dynamo action is observed, is small. 
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1. Introduction 



The present work continues the series of studies of generation of magnetic field pos- 
sessing large spatial scales by Lanotte et al. (2000), Zheligovsky et al. (2001) and 
Zheligovsky and Podvigina (2003). Space-periodic parity-invariant flows, steady or 
periodic in time, which have an exponentially decaying spectrum, were considered by 
these authors, regarded as a model for natural turbulent flows of conducting fluids. 
For such flows magnetic eddy diffusivity (MED) was found to be often negative. (The 
analysis has been extended to the study of linear stability of space periodic magne- 
tohydro dynamic steady states to long-period perturbations by Zheligovsky, 2003). 

We consider here a kinematic dynamo problem in an infinite layer, assuming 
perfectly conducting horizontal boundaries: 2 



dH l 



213 = 0, 7T 



= 0, #3|s3=0,7T = 0. (1) 

X3=0,7r 



As ibid., the following algebraic idea serves as a foundation for our constructions. 

The kinematic dynamo problem for a steady flow v(x) can be reduced to an 
eigenvalue problem for the magnetic induction operator C: 

CH = r]V 2 H + V x (v x H) = AH. (2) 

The adjoint operator is 

£*H = 7]V 2 H - v x (VxH), 

implying £*H = for any constant vector field H. Hence, any constant vector field, 
satisfying the boundary conditions for the adjoint operator, belongs to its kernel. (By 
the definition of the adjoint operator, the identity 

(£H 1 ,H 2 ) = (H 1 ,£*H 2 ) (3) 

holds for all vector fields Hi(x) and H 2 (x) from domains of C and £*, respectively; 
here the standard scalar product (•, •) of the functional Hilbert space L 2 (the Lebesgue 
space) is assumed. Boundary conditions for the adjoint operator can be determined 
demanding that all surface integrals appearing in (3) after integration by parts van- 
ish.) Generically the kernel of C* is spanned by such constant vector fields. Thus 
generically there exist as many neutral magnetic modes (i.e. magnetic fields, sat- 
isfying (2) with A = 0), as there exist linearly independent constant vector fields 
satisfying the boundary conditions for C*. 

The flow is supposed to be periodic in horizontal directions: 

V(x) = V (xi + jr-,X2, X^j =V (x 1 ,X 2 + X 3 ^j . (4) 

If vector fields in the domain of C have the same periodicity (4) as the flow, 

2tt \ „ / 2tt 



H x 



H (xx + -j^-,x 2 ,x 3 ^ = H (x u x 2 + J^i x z) » ( 5 ) 



2 Roberts and Zhang (2000) commented on physical validity of these boundary conditions. 



2 



(vector fields with such periodicity will be called small-scale) and satisfy boundary 
conditions (1), then vector fields in the domain of C* also satisfy (5) and (1), provided 
the fluid does not penetrate through horizontal boundaries of the layer: 

V 3 \x 3 =0,n = 0. (6) 

Hence in this case a constant vector field belongs to the kernel of C*, if and only if 
its vertical component vanishes. Thus, 

dim ker C > 2 

and generically there exist two neutral magnetic modes satisfying (5) and (1). 

The conducting fluid resides in an infinite volume. When the spatial period of 
magnetic field is allowed to increase to infinity, the smallest eigenvalue of the Lapla- 
cian tends to zero - diffusive dissipation of slowly varying magnetic fields is small. 
Therefore, the following question is legitimate: Can a growing magnetic mode be 
constructed by perturbing a neutral mode of the same periodicity as that of the flow, 
allowing larger spatial scales in the perturbation? It was found in the studies cited 
above, that growing modes can indeed be constructed following this idea and us- 
ing multiscale techniques. We show in the present paper, that this remains true for 
magnetic field generation in a layer. 

The mechanism of negative MED for generation of large-scale magnetic field re- 
quires that no a-effect is present. To ensure this, parity-invariant flows have been 
considered in the earlier studies. In the present work a different symmetry is employed 
for this purpose - the symmetry about the vertical axis: 

v 1 (-x 1 ,-x 2 ,x 3 ) = -v 1 (x 1 ,x 2 ,x 3 ), 

v 2 (-x 1 , ~x 2 , x 3 ) = -v 2 (x 1 ,x 2 , x 3 ), (7) 

V3(~Xl, -X 2 ,X 3 ) = V 3 (X!,X 2 ,X 3 ). 

This type of symmetry is not unnatural: it is sustained by the Navier-Stokes equation 
(if the external forces involved have this symmetry). 

In Section 2 we discuss expansions of magnetic modes and their growth rates 
in power series in scale ratio, e, which is a small parameter of the problem, and 
present MED tensor for a general incompressible flow in a layer satisfying (4) and 
(6). A complete formal expansion of magnetic modes and the associated eigenval- 
ues is exposed in detail in Appendix A. In Section 3 we prove that for plane flows 
(and hence for convective rolls) eddy correction to molecular magnetic diffusivity is 
always non-negative. In the remaining part of the paper results of Section 2 are ap- 
plied to convective plan-forms without rotation. Relevance of the kinematic dynamo 
problem for convective plan-forms for the study of the full non-linear magnetohydro- 
dynamic system has been discussed in the context of the theory of bifurcations by 
Bosh-Vivancos 3 et al. (1995). Representation and symmetry properties of convective 

3 Numerical results of this paper have been questioned by Zheligovsky and Galloway (1998) and 
Matthews (1999b). 
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plan-forms are summarised in Section 4. The plan-forms possess additional symme- 
tries, for which the MED tensor is diagonal, as discussed in Section 5. In Section 6 
numerical results for two families of plan-forms are presented. For cell patterns pos- 
sessing the symmetries of a rectangle, the region in the parameter space of negative 
MED coincides with the region of small-scale magnetic field generation. A family of 
cells with a smaller symmetry group has been found numerically, where MED is neg- 
ative for molecular magnetic diffusivity over the threshold for the onset of small-scale 
magnetic field generation. In Appendix B we consider alternative representations 
of elements of the MED tensor and procedures for computation of MED, which are 
more efficient than straightforward numerical solution of auxiliary problems arising 
in construction of the asymptotic series. 



2. Magnetic eddy diffusivity tensor 

for large-scale magnetic modes in a layer 

We outline here construction of expansions of magnetic modes and their growth rates. 
This is done under the basic assumption that a generic problem is considered, i.e. the 
kernel of the adjoint operator C* is spanned by constant vector fields, whose vertical 
component vanishes. In this Section a steady fluid flow v(x) is supposed 

i. to satisfy the boundary conditions (4) and (6); 

ii. to possess the symmetry (7); 

Hi. to have a zero space average of the horizontal components; 
iv. to be solenoidal: 

V • v = 0. (8) 

A magnetic mode H(x, y) is supposed to depend on the fast variable x e R 3 and 
on the slow variable y = e(rri,£ 2 ) in horizontal directions. Here the scale ratio e > 
is a small parameter. A solution to the eigenvalue problem (2) is sought in the form 
of power series in e: 

oo 

H=5> (n) (x,y)e», (9) 

n=0 

oo 

A = A (ri) e ra . (10) 

?1=0 

The mode is supposed to be 27r/L. ; -periodic in Xj for j = 1,2, and to satisfy (1) on 
the horizontal boundaries x% = and x 3 = it. Accordingly, we assume that each h^) 
satisfies (5) and (1). 

Let (•) and (( • )) denote the mean part of a scalar or vector field and the mean 
horizontal part of a vector field, respectively: 




Here V = 47r 3 /(L 1 L 2 ) is the volume of the periodicity box, and ej is a unit vector in 
the direction along the coordinate axis Xj. Denote 

H (n) = ^ h (n)y 
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By the chain rule, dependence of the magnetic mode on the fast and slow variables 
implies that gradients must be modified in the eigenvalue equation (2): 



V x + eV y , where V y = 



— — o" 

dyi ' dy 2 ' 



(11) 



Substituting (9)-(ll) into (2), one obtains a hierarchy of equations, which can be 
solved successively in all orders of e (see Appendix A). 

In particular, one finds A = 0(e 2 ) : A^ = \^ = 0, and A^ is determined 
from the eigenvalue equation for the mean horizontal part of the leading term in the 
expansion of the magnetic mode (9): 

2 2 orr(O) 

MH^ = 77VjH(°) +PV y x^^(vx r ( m , k)) ^k_ = A (2) H (o)_ (12) 

k=i m=i oy m 

Here V is the projector onto horizontal directions: 

j'=i 

Ai is called the operator of (anisotropic) MED. Its coefficients can be determined 
from two auxiliary problems: 

the first auxiliary problem: 

£S (fc) = -f^ (A; = 1,2); (13) 
the second auxiliary problem: 

CT (m,k) = _ 2 e m x (v x (SW + e fc )) (m, A: = 1,2). (14) 

Vector fields S^ fe ^ and r( m,A: ^ satisfy the boundary conditions (5) and (1). The basic 
assumption (stated in the beginning of this Section) implies solvability of the problems 
(13) and (14). It is shown in Appendix A that 

V-S (fc) = 0, (15) 

V • T {m ' k) + = 0. (16) 

Alternative expressions for elements of the MED tensor can be obtained, partially 
performing integration required for averaging of the cross-products of v and T^ m ' k \ 
They are derived in Appendix B. 

M. is a second order operator in partial derivatives with constant coefficients. 
Consequently, solutions to (12) bounded on the entire plane are Fourier harmonics 

H (0) = he iq-y ; 
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where q G R 2 is an arbitrary constant wave vector, and h G i? 3 is a constant vector 
satisfying 

77|q| 2 h + V ([q u q 2 , 0]x^(vx T^)q m h k ) = -A^h (18) 

V fc=lm=l / 

and 

h 3 = 0. (19) 
The solenoidality condition for the magnetic mode 

V • H = (20) 

implies V y • H^ ^ = (see Appendix A), hence 

[ qi ,q 2 ,0}-h = 0. (21) 

From (18), (19) and (21), 



2 2 



h= [-92,91,0], -A< 2 >(q) =vh\ 2 -Y: E(-lf(vx T^%q m q^ k . (22) 

k=l m=l 

If |q| = 1, following the long established tradition — is called magnetic eddy 
diffusivity. If the minimal MED is negative: 

min(-A( 2 )(q)) < 0, 

|q|=i 

there exist growing large-scale magnetic modes, i.e. the flow operates as a dynamo. 
The leading term of the magnetic mode expansion (9) can be expressed as 



h(°)(x,y)=e i ^5:/ ifc (S( fc )(x)+e fc ). 



k=i 



Since (13) is equivalent to 

£(SW + e fc ) = 0, (23) 
the magnetic mode is a perturbation of the neutral small-scale magnetic mode 

f> fc (S (fe) (x) + e fc ) Gk er£ 
k=i 

modulated by an amplitude factor depending on the slow variable. 



3. Magnetic eddy diffusivity for plane parallel flows 

By the Zeldovich (1956) theorem (see also Moffatt, 1978), plane flows (such flows 
v(x) that v(x) ■ k = for a constant vector k) cannot generate magnetic field. Here 
a stronger in some sense result is derived: plane parallel flows can only enhance 
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molecular diffusivity. More precisely, we show that magnetic eddy correction for plane 
parallel flows in a layer, which satisfy conditions i — iv stated in the first paragraph 
of Section 2, is always non-negative. Without any loss of generality we assume in this 
Section, that the flow does not depend on x 2 and v 2 = 0. 
k = 1. From (13), 

Si 1] = v x (S (1) + ei) = i?(x)e 2 . 

Thus for m — 1 the right-hand side of (14) equals —2rfdS^/dxi — i?(x)e 3 ; conse- 
quently 

rf' 1} = o (v x r* 1 - 1 ^ = o. 

For m = 2, the right-hand side of (14) vanishes and hence 

r^ 2 ' 1 ) = o (v x r^ 1 ^ = o. 

k = 2. From (13), = 0. Hence, for m — 1 the right-hand side of (14) is i>ie 2 . 
This implies I^ 1 ' 2 ) = re 2 , where T satisfies 

?7V 2 r- (v v)r = «i. 

A similar equation arises in the multiscale analysis of the passive scalar transport 
equation (see Biferale et al. , 1995). The respective coefficient of the MED tensor is 
non-negative: 

-(vxr( 1 ' 2 )) 3 = -(t; 1 r) = ^(|vr| 2 )>o 

(it is positive for v^O). 

For m — 2, the right-hand side of (14) is equal to — v, and hence 

if ' 2) = (v x r( 2 < 2 )> 3 = o. 

Consequently, (22) implies 

-A( 2 )=77|q| 2 -(vxr( 1 ' 2 )) 3g 2 >?7|q| 2 , 

and thus minimal MED is equal to molecular magnetic diffusivity rj, as we intended 
to demonstrate. 

4. Convective plan-form flows in a layer 

Analysis of Section 2 will be applied to convective plan-forms without rotation. In this 
Section we summarise some of their properties (see Chandrasekhar, 1981). Thermal 
convection in a fluid heated from below in a layer with no rotation is considered. 
Plan-forms are instability modes (more precisely, the flow parts of the modes) at the 
onset of instability of the trivial steady state (in which the fluid is at rest and the 
temperature profile is linear in x 3 ). They are poloidal: 

v(x) = V x V x (P(x)e 3 ). (24) 
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In the most general form (Bisshopp, 1960) the potential can be expressed as 

P(x) = (cti cos(LiXi) cos(L 2 ^2) + a 2 cos(pL 2 x 2 ))w(x 3 ). 
Here a\ and a 2 are constant and p is integer. For a 2 ^ 0, 



(25) 



uiu = - 1. 



Evidently, plan-forms satisfy conditions ii — iv assumed in Section 2. Fluid does 
not penetrate through the horizontal boundaries (condition i), if and only if 



w(x 3 )\ 



X3=0,TT 



= 0. 



(26) 



If in addition 



d 2m w 



dx 2m 



Vm > 0, 



the plan-form satisfies the free boundary conditions: 



dvo 



X 3 =0,7T 



0, U3Us=0,ir = 0- 



X 3 =0,TT 



Consequently, for free boundaries the vertical profile is 

w n (x 3 ) = sinnx 3 , 



(27) 



where n > is integer. The mode for n — 1 is the first to become unstable when the 
Rayleigh number is increased. For rigid horizontal boundaries the vertical profile is 
more complex (see Chandrasekhar, 1981). 

Plan-forms possess two "reflection" symmetries, referred to in the next Section: 
reflection in the direction ei: 



(28) 



l>l(-Xl, 


x 2 ,x 3 ) = 


-v 1 (x 1 ,x 2 ,x 3 ), 


V 2 (~Xi 


, x 2 , x 3 ) = 


v 2 (x 1 ,x 2 ,x 3 ), 


V 3 (-X! 


,x 2 ,x 3 ) = 


v 3 (x 1 ,x 2 ,x 3 ); 


e 2 : 






vi(x u - 


-x 2 ,x 3 ) = 


v 1 (x 1 ,x 2 ,x 3 ), 


v 2 (xi, - 


x 2 ,x 3 ) = 


-v 2 (x 1 ,x 2 ,x 3 ), 


v 3 (x 1: - 


-x 2 ,x 3 ) = 


v 3 (x u x 2 ,x 3 ). 



(29) 



(The symmetry about the vertical axis is a composition of these two.) 

Free boundary plan-forms with the vertical profile (27) are parity-invariant about 
centers specified in the following table: 
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odd 



odd 



even 



even 



a 2 = 



odd, or a 2 = 



any 



even, or «o = 



71 



Centers 



nl 2 7r N 
'V2, 



I Till 7ll 2 7r" 

1 17' V 2, 



7T 

' 2 



h + 



1 



7T 



Here l\ and Z2 are arbitrary integers. For such plan-forms, the same asymptotic 
expansions could be constructed using parity invariance instead of (7). 



5. The operator of magnetic eddy diffusivity 
for convective plan- form flows in a layer 

Each of the two reflection symmetries (28) and (29) splits the domain of £ into a 
direct sum of two proper subspaces: symmetric or antisymmetric with respect to the 
given symmetry. (A vector field is said to possess the reflection antisymmetry in the 
direction ei, if it satisfies (28) with the reversed signs in the right-hand side: 

v 1 (-x 1 ,x 2 ,x 3 ) = Vi(xi,X 2 ,X 3 ), 

V 2 (-X!,X 2 ,X 3 ) = -v 2 (x!,x 2 ,x 3 ), 
v 3 (-x 1 ,x 2 ,x 3 ) = -v 3 (x 1 ,x 2 ,x 3 ); 

a reflection antisymmetry in the direction e 2 is defined similarly.) 

Consequently, due to (13) the solution to the first auxiliary problem has the 
same reflection symmetries/antisymmetries, as the vector field dv/dxk- Furthermore, 
due to (14) the solution r*™'^ to the second auxiliary problem has the same reflection 
symmetries/antisymmetries, as d 2 v/dx m dxk. 

In particular, r( fc ' fc ) possess both reflection symmetries and hence 

(v x r (fe ' fc) ) 3 = 0. 

This implies diagonality of the MED operator, restricted to the subspace of solenoidal 
vector fields of slow variables with the zero vertical component: 

M = ( V - (v x r^> 3 )|^ + (v + (v x r^s) jg. 

Therefore, for a flow possessing the symmetries (28) and (29), the minimal MED is 

min(-A (2) (q)) = min^ - (v x T (1 ' 2) ) 3 , r] + (v x r {2 ' 1} ) 3 ). (30) 
|q|=i 
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6. Numerical results 



As shown in the previous Section, to find coefficients of the MED tensor and thus the 
minimal MED for a plan-form it appears necessary to solve four problems £X = F: 
two first, and two second auxiliary problems (13) and (14). However, as discussed in 
Appendix B, it suffices to solve three elliptic problems of the same complexity. 
Computations have been performed for flows (24), (25) with the vertical profile 

N 

w ( x 3) = £ A. sin(nx 3 ), (31) 

n=l 

where f3 n are constant coefficients. The flows employed in computations are nor- 
malised, so that the root mean square is equal to 1. All Figures below show results 
obtained for plan-forms with the vertical profile (27). 



6.1 Fourier representation of the vector fields and r^ m fc ^ 

Since solutions to the auxiliary problems possess the reflection symmetries (28) and 
(29) (see Section 5) and satisfy the boundary conditions (1), they can be expanded 
in the following Fourier series: 



S (D 



E 

n;>0 



s ni cosniLi^i cosn 2 L 2 x 2 cos 713X3 

(i) • T T 

s n2 smn i^i x i smn 2 L 2 x 2 cosn 3 x 3 

(!) T T 

s„ 3 smniLiXi cosn 2 L 2 x 2 sin 713X3 



E 

n;>0 



s (2) 
,(2) 



smriiLiXi sinn 2 L 2 x 2 cosn 3 x 3 
Sn2 cosn i^i x i cos n 2 L 2 x 2 cosn 3 x 3 

(2) T ■ T 

Sn 3 cosniLiXi sin n 2 L 2 x 2 sin 713X3 

(m,k) 



(32) 



(33) 



7nT costiiLiXi sin n 2 L 2 x 2 cos 713X3 
p(m,fc) _ ^2 sinriiLiXi cosn 2 L 2 x 2 cosn 3 x 3 . (34) 

n;>0 (m,k) . T -j 

L 7n,3 smniLiXi sm n 2 L 2 x 2 sm 713X3 

All computations have been performed with the resolution of 64 trigonometric 
functions in each direction. With this resolution the energy spectrum of solutions 
decays by several (at least 4) orders of magnitude. 

The following parity symmetries each reduce twice the number of unknown coef- 
ficients. If in (25) p is odd or a 2 = 0, the sum of horizontal wave numbers in the 
velocity is even; hence vector fields constructed of harmonics, where the sum of the 
wave numbers in the horizontal directions is either even or odd, constitute, respec- 
tively, two invariant subspaces of the domain of C Therefore, if p is odd or a 2 = 0, 
coefficients of the series (32)-(34) for odd n\ + n 2 vanish. Similarly, if in (31) only 
odd wave-number terms are present, coefficients of the series (32)-(34) for odd rii + n 3 
vanish. 

For plan- forms with the vertical profile (27) for n > 1, coefficients of the series (32)- 
(34) vanish for all n 3 , which are not divisible by n, and the second parity symmetry 
is modified: coefficients are non-zero only for even ri\ + n 3 /n. 



10 
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Figure 1. Minimal MED (solid line; vertical axis) and the growth rate of the dominant 
small-scale magnetic mode with a zero horizontal mean (dashed line; vertical axis) 
for cell patterns possessing the symmetries of a rectangle, for L 2 = 2, rj = 0.06, as a 
function of L x (horizontal axis). Computed values are shown by solid dots. 

6.2 MED for rectangular cell patterns 

Thermal convection cell patterns possessing the symmetries of a rectangle are plan- 
forms (24), (25), where a 2 = 0. 

Matthews (1999a,b) 4 found square plan- forms with L\ — L 2 — 1/2 incapable of 
kinematic dynamo action. He assumed the boundary conditions (1) and considered 

4 Opposite to Matthews (1999b), we cannot recommend integration of the magnetic induction 
equation in time as an efficient technique for numerical determination of dominant magnetic modes 
for steady flows, at least unless the following improvements are implemented: 1) Integration of 
the magnetic induction equation in time is performed using fast specialised time-stepping schemes 
overcoming stiffness of the problem (e.g., see Nikitin 1994, 1996). (They can be easily implemented, 
if the Galerkin discretisation of the induction equation in space in the basis of eigenfunctions of the 
Laplacian is employed.) 2) Optimisation methods (e.g., similar to those of Zheligovsky, 1993) are 
employed, allowing one to jump from the current trajectory of temporal evolution of a magnetic 
field to a different one, closer to the exponential trajectory for the dominant magnetic mode. 3) 
Smallness of the discrepancy |£h — £h|, and not the "overall exponential behaviour" of the obtained 
solution is examined in the condition for termination of computations. 

In the present work a spectral code derived from that of Zheligovsky (1993) has been used. 
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Figure 2. Minimal MED (solid line; vertical axis) and the growth rate of the dominant 
small-scale magnetic mode with a zero horizontal mean (dashed line; vertical axis) 
for cell patterns possessing the symmetries of a rectangle, for Li = 0.75, L 2 = 2, as 
a function of rj (horizontal axis). Computed values are shown by solid dots. 

magnetic modes with vectors of basic periods along the diagonals of our basic peri- 
odicity square, with the periods equal to a half of the diagonal length (constituting a 
subspace of our small-scale magnetic fields). 

Figures 1 and 2 display MED computed for two sets of parameter values for plan- 
forms with the symmetry group of a rectangle, together with the growth rate of the 
dominant magnetic mode of the same spatial periodicity as in the flow (the associated 
eigenvalues of C are real). 

The Figures illustrate an unexpected phenomenon: The critical molecular viscos- 
ity for the onset of small-scale magnetic field generation coincides with that for the 
onset of large-scale magnetic field generation; i.e. the minimal MED vanishes together 
with the growth rate of the dominant magnetic mode of the spatial periodicity of the 
flow. Since the growth rate of the large-scale magnetic mode is small (0(e 2 )), the 
large-scale magnetic instability is weak compared to the small-scale one. 

This phenomenon has been reproduced in simulations for a number of vertical 
profiles (31) comprised of up to three sines with wave numbers of the same or different 
parity. In these simulations the critical molecular diffusivities r) have been determined 
to the accuracy of at least 10~ 8 . Algebraic reasons of this phenomenon are not entirely 
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Figure 3. Minimal MED (solid line; vertical axis) and the growth rate of the dominant 
small-scale magnetic mode with a zero horizontal mean (dashed line; vertical axis) 
for plan-forms for p = 2, U = 1 V3, L 2 = 1.5, 7] = 0.0075 as a function of a 2 /a\ 
(horizontal axis). Computed values are shown by solid dots. 



clear. 

Suppose horizontal coordinate axes are such that 77 — (v x I"^ 1,2 )^ is the minimal 
MED (cf. (30)); they can be always chosen this way, since for a 2 = the two horizontal 
coordinate directions are interchangeable. Then the neutral small-scale magnetic 
mode has the same reflection symmetries as I^ 1 ' 2 ), and can be expanded in a Fourier 
series of the form of (34). All non- vanishing terms in this series correspond to wave 
vectors, whose n components have the parity either (odd, even, even), or (even, odd, 
odd). This is in contrast with the structure of I^ 1,2 ), which is represented by the series 
(34) where all non- vanishing terms correspond to wave vectors, all three components 
of which have the same parity. 

6.3 MED for general cell patterns 

Figure 3 displays MED computed for a family of cell patterns for a 2 7^ (which do 
not possess therefore the symmetries of a rectangle) together with the growth rate 
of the dominant magnetic mode of the same spatial periodicity as in the flow (the 
associated eigenvalues of £ are again real). A window of small ratios a 2 /a± has been 
detected, where MED is negative. For the assumed value of molecular viscosity flows 
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from this family do not generate small-scale magnetic field; thus in this family of 
plan-forms the effect of negative MED instability is separated out. 

For «2 = cki/2 the flow is the Christopherson (1940) hexagonal cell pattern. Zhe- 
ligovsky and Galloway (1998) found that in case there is a dielectric beyond one of 
the horizontal boundaries of the layer, the flow can generate small-scale magnetic 
field for small molecular diffusivities. Matthews (1999b) has not found magnetic field 
generation by hexagonal cells for L x = \j3/8, L 2 = \Jl/&, assuming the boundary 
conditions (1) and (5). A number of computations for the hexagonal cell pattern have 
been performed for two-term sums (31) and a varying ratio f3 n / Pi, but no instances 
of negative MED have been found. 



7. Conclusion 

Formal asymptotic expansions of magnetic modes and the associated eigenvalues in 
series in the scale ratio have been constructed for flows, symmetric about the vertical 
axis. The region where the flows reside and their symmetries are different from those 
for parity-invariant space-periodic flows, which have been considered so far. If a flow 
possesses reflection symmetries about the horizontal Cartesian axes, the magnetic 
eddy diffusivity tensor becomes diagonal. The suggested strategies for evaluation 
of elements of the MED tensor require numerical solution of a smaller number of 
elliptic PDE's, than for the direct approach whereby all auxiliary problems are solved. 
Numerical experiments have been performed for convective plan-forms in the absence 
of rotation, and a family of plan-forms exhibiting negative MED has been found; 
these flows cannot generate small-scale magnetic field with the spatial periodicity of 
the flow. Therefore magnetic field generation in a layer by the mechanism of negative 
MED is possible. The flows are poloidal, and their helicity vanishes identically. Thus, 
the paradigm requiring helicity of a flow to be non-zero to enable it to act as a dynamo 
has been once again demonstrated to be wrong. However, the region in the parameter 
space, where large-scale dynamo operates, is small, and the minimal negative MED's 
found in computations are small in absolute value. 
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A. Asymptotic expansion of magnetic modes 
and the associated eigenvalues for flows, 
symmetric about the vertical axis 

Complete formal asymptotic expansions of magnetic eigenmodes and associated eigen- 
values are constructed here in the form (9), (10) for small values of the scale ratio 
e. Our basic assumption is that the problem is generic, i.e. the kernel of the adjoint 
operator C* is spanned by constant vector fields, whose vertical component vanishes. 
The remaining assumptions concerning the flow are listed in the first paragraph of 
Section 2. 

Evidently, ((£X)) = 0. Let us prove that the basic assumption is equivalent 
to the assumption that C, restricted to the subspace of vector fields which have a 
zero horizontal mean and satisfy (5) and (1), is invertible. An equation £X = F is 
solvable, if and only if F is orthogonal to vector fields from ker C*. By the assumption, 
if ((F)) = 0, this condition is satisfied. Let us show that X satisfying ((X)) = can be 
found. Suppose there exists a e ker£, ((a)) = 0. C* is an elliptic operator; for the 
boundary conditions (5) and (1) it has a point spectrum and its eigenvectors a£, 

= £^0forfc>3, 

constitute a complete basis (assuming for the sake of simplicity that it does not have 
generalised eigenvectors). Expand a in this basis: 

E* \ ^ Cfe p* * 

k>3 fc>3 ^fc 

Scalar multiplying this equation by a one finds |a| 2 = 0. (Alternatively, an argument 
similar to that of Jones and Roberts, 2000, can be put forward: averaging a horizontal 
component of the equation for a small-scale magnetic mode, 

for the assumed boundary conditions one finds = £fc(( a fc)) 5 implying that all small- 
scale magnetic modes with a non-zero horizontal mean belong to ker£.) Thus any 
non-zero vector field from the kernel of C has a non-zero horizontal average. Conse- 
quently, it is possible to enforce ((X)) = 0, subtracting from X the appropriate vector 
from ker C. This concludes the demonstration of equivalence of the two assumptions. 
Denote (in accordance with Section 2) 

H W = (( h (n)^. Q(n) = h (n) _ ^(n)^ 

A mode is supposed to be L 3 -periodic in Xj for j — 1,2 and to satisfy (1) on the 
horizontal boundaries. Consequently, we assume that each h( ra ) in (9) satisfies (5) 
and (1). Since ((h^)) automatically satisfies them, G^ also does. 

Substituting the series (9) and the expression for the spatial gradient (11) into 
the solenoidality condition (20) and expanding the result, find 

oo 

(V y • H (n - X) + V x • G (n) + V y • G (n_1) ) e n = 

n=0 
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(by definition, H^ n ^ = = Vn < 0). This defines a hierarchy of equations. The 
mean of the equation at order n + 1 is 

V y • H("» = 0; 

the equation at order n after subtraction of the mean becomes 

V x • G (n) + V y • G (n ' 1} = 0. 

Substituting (9)-(ll) into the eigenvalue equation (2), obtain its series represen- 
tation: 



E 

n=0 



£G (n) + r] (2(V X • V y )G (n " 1) + V y (H( n " 2 ) + G ( "- 2) )) 
+V X x (v x H (n) ) + V y x (v x (H (n_1) + G ( "- 1} )) (A.1) 



e n = 0. 



— V] \( n ~ m )(ii( m ) _|_ Q( m )^ 

m=0 

This yields a hierarchy of equations of the form 

£G (n) = F (ri) (x,y). (A.2) 

By the basic assumption this equation is solvable for G^ if and only if ((F^ (x, y))) = 0. 

Equations (A.2) arising from (A.l) can be solved successively in all orders by the 
following two-step procedure: 

1. for a given n, from the solvability condition H^™ -2 ) is determined and the solvability 
condition is thus satisfied; 

2. the resulting equation for this n is solved for G^ n \ 
i. n = 0. (A.2) reduces to 

CG (o) + ( H (o) . Vx)v = A (o) (H (o) + G ( 0)) _ 

The mean horizontal part of this equation is = A^H^) implying = (we 
consider modes which are not predominantly small-scale, i.e. such that H^ ^ ^ 0). 

By linearity, the solution with a zero horizontal mean to the resulting equation 
admits a representation 

G(°) = £ S ( fc )(x)4 0) (y), (A.3) 

k=l 

where S^(x) are solutions to the first auxiliary problem (13), satisfying (5) and (1). 
The solvability condition for (13) is verified, since { — dv/dxk} = due to periodicity 
of v. In view of (23) equivalent to (13), + e k span the kernel of C 

Since v is symmetric about the vertical axis, the domain of £ is a direct sum of 
two proper subspaces, one of which is comprised of vector fields, symmetric about the 
vertical axis: 

H 1 (-x 1 ,-x 2 ,x 3 ) = -H 1 (xi,x 2 ,x 3 ), 
H 2 (-x 1 ,-x 2 ,x 3 ) = -H 2 (x 1 ,x 2 ,x 3 ), 
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H 3 (-x ± , -x 2 , x 3 ) = H 3 (x!,x 2 , x 3 ); 
and the second - of the ones, antisymmetric about the vertical axis: 

H 1 (-x 1 , -x 2 ,x 3 ) = H 1 (x 1 ,x 2 ,x 3 ), 

H 2 (-X!, -x 2 ,x 3 ) = H 2 (X!,X 2 ,X 3 ), 

H 3 (-x 1 ,-x 2 ,x 3 ) = -H 3 (x!,x 2 ,x 3 ). 

Consequently, is antisymmetric about the vertical axis. 
Divergence of (13), 

V 2 (V- S (fc) ) = 0, (A.4) 

suggests solenoidality of S^ k \ Because of boundary conditions (1) for and (6) for 
the flow, the vertical component of (13) implies d 2 S 3 k ^ /d 2 x 3 = 0. Consequently, 

£3 = 0, 7T 

V • SW satisfies 

^-v • s( fc ) = 

ox 3 

on horizontal boundaries. Evidently it has the same periodicity (5) as the flow. 
Multiplying (A.4) by V ■ and integrating over the periodicity box, one obtains 
(|V(V • S (fc) )| 2 ) = =>- V ■ S (/c) = constant. Integration of this equation over the 
periodicity box demonstrates solenoidality of S^ h \ 
ii. n = 1. (A. 2) takes now the form 

£G« + 2r?(V x • V y )G<°> + (H« • V x )v 

+V y x (v x G(°)) - (v • V y )H(°) = A( 1 )(H(°) + G (0) ). 
The mean horizontal part of this equation is 

PVyX^vxS^f'^WHf '. 
k=i 

Since v is symmetric about the vertical axis and S^ fc ^ is antisymmetric, 

(v x S (fc) ) 3 = A (1) = 0. 

Thus, although in the system under considerations the a-effect tensor does not vanish 
entirely, it does not affect generation of the principal component of the mean field. 
Consequently, after (A. 3) is substituted, the equation becomes 



2 2 / 

£G (D = _ (H « ■ Vx )v -EE 2t^- + e m x (v x (S<*> + e fc )) 

fc=lm=l V ° Xm / 



dy m ' 

By linearity, its solution with a zero horizontal mean admits a representation 

G (i) = E S <*>(x)ltf>(y) + E E r^(x)^-(y), (A.5) 

k=i fe=im=i °y m 
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where vector fields T^™ 1 '^ are solutions to the second auxiliary problem (14) satisfying 
the boundary conditions (5) and (1). Since v and are symmetric and antisym- 
metric about the vertical axis, respectively, and ((v)) = 0, the average of horizontal 
components of the right-hand side of (14) vanishes, implying solvability of (14). 

Due to the antisymmetry of and the symmetry of v, the right-hand side of 
(14) is symmetric about the vertical axis, and since £ preserves the symmetry and 
antisymmetry, r*" 1 '^ are symmetric about the vertical axis. 

Let us show (16). Denote by $ a vector potential of v: 

v = V x 



From (13), 

v x S (/c) = ??V x S (/c) — — h V^ (fc) + C (k) . (A.6) 

Here <1> and ip^ have the flow's periodicity (4), and is a constant vector. Now 
(14) can be represented as 

cr(m,k) + v [~dx~ + vs ™) + v x (($fe " ^ {k))ern) + GmX c(k) = °- 

Taking divergence of this equation, find 

V 2 (V • T (m ' fc) + S^) = 0. (A.7) 

Because of boundary conditions (1) for and Y^ m ' k \ and (6) for the flow, the ver- 
tical component of (14) implies d 2 T ( £ n ' k ' 1 /d 2 X3 = 0. Consequently, the quantity 

.i':; II.:: 

= V • r*™'^ + Sj® satisfies dcft/dx^, = on horizontal boundaries. Evidently, it has 
the same periodicity (5) as the flow. Multiplying (A.7) by <fi and integrating over the 
periodicity box, obtain (|V0| 2 ) = =>- <fi = constant. Integration of this equation 
over the periodicity box demonstrates (16). 
in. n = 2. (A. 2) takes the form 

CG™ + V (2(V X • V y )G« + V 2 (H<°> + G(°))) 

+ (H< 2 > • V x )v + V y x (v x G«) - (v • V y )H« = A^tK ) + G (0) ). 

The mean horizontal part of this equation is (12). Its solution (17), (22) has been 
derived in Section 2. Thus G^ is determined together with H(°) from (A. 3). The 
equation for G^ can now be solved by the general procedure outlined in iv for n > 2. 
iv. n = N > 2. From equations for n < N one finds: 

• vector fields G^ n ^ for n < N — 2, in particular 

GW(xj)=g( n '(x)e iw ; (A.8) 

• H (n) = for < n < N - 2; (A.9) 
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• representations for of the form 

G W = J- S«(x)^ n) (y) + E E r< m ' fc >(x)^ (y) + Q^(x)e iq ' y (A. 10) 

k=l k=l m=l ^ m 

for n = iV - 1 and n = N - 2 with known vector fields Q (n) , ((Q (ri) )) = (note (A.3) 
and (A. 5) are particular cases of (A. 10) with Q(°) = = 0); 

• A< n > for n < N. 

The mean horizontal part of (A. 2) is 

N-3 

(M - A (2) )H^- 2 ) - A (JV) H(°) = -V{V y x (v x Cp- 1 ))) + ]T \ {N ~ m) H^; (A.ll) 

m=l 

its right-hand side is known. In view of (A.8)-(A.10), dependence of the right-hand 
side of (A.ll) on the slow variable is via the factor e iqy . Projecting this equation 
out in the direction of H( ) in the space of bounded solenoidal vector fields of slow 
variables, whose vertical component vanishes, one uniquely determines A^. 

In the complementary invariant subspace in the space of solenoidal vector fields 
of slow variables only, whose vertical component vanishes, the operator M. — A*- 2 - 1 
is invertible. Consequently, H^ -2 ) = up to an arbitrary multiple of H^°\ which 
can be neglected (this is a normalisation condition). Therefore, G^ -2 -* is now also 
determined. 

The resulting equation for G^ becomes 

£G<"> = -r, (2(V X ■ Vj-JQ""- 11 + VjG""- 21 




dx m dy rn fr{ dxi dyidy. 

-v y x ( E(v x ^)Hi N ' l) + E E { v x r(m ' fc) } ^7— + { v >< Q (iV ~ 1} } 



gjjj.N-2) 



v fe=l k=l m=l 

N-2 



+ E A (Ar ~ m) G (m) , 

rn=0 

where it is denoted {f} = / — ((/)). Only terms involving and derivatives of 

Hl N ~^ have not yet been determined in it. x-dependent pref actors in front of these 
terms are identical to those in front of and derivatives of , respectively, in the 
equations for n = and n — 1. Therefore, G^ admits the required representation 
(A. 10). (An equation for is obtained by replacing by and omitting 
all terms, involving and derivatives of H^ -1 ).) 

Furthermore, it can be shown that Q^* 1 and G^ are symmetric about the vertical 
axis for odd N, and antisymmetric for even N. Consequently, from (A.ll), 

V(V y x (v x Q^^)) = \ {N) = for odd N. 

Thus complete formal expansions of magnetic modes (9) and the associated eigen- 
values (10) have been constructed; they satisfy (17), (22), (A. 9) and (A. 8). 
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B. Computation of elements 

of the magnetic eddy diffusivity tensor 

Here we discuss strategies for evaluation of elements of the MED tensor, requiring 
to solve a smaller number of elliptic PDE's than for the direct approach whereby all 
auxiliary problems are solved, and present alternative expressions for elements of the 
MED tensor, where integration required for averaging of the cross-products of v and 
r(m,fc) j s partially performed. 

As derived in Appendix A, elements of the MED tensor are equal to (v x T^ m ' k ^)i. 
For flows symmetric about the vertical axis in a layer, considered in this paper, I = 3 
and m, k = 1,2. For parity-invariant space-periodic flows, auxiliary problems are 
represented by the same equations (13) and (14), but now l,m,k = 1, 2, 3. To evaluate 
the MED tensor, the first and second auxiliary problems were solved by Lanotte et 
al. (2000), Zheligovsky et al. (2001) and Zheligovsky and Podvigina (2003). This 
approach requires solution of 6 elliptic PDE's £X = F for flows symmetric about the 
vertical axis in a layer, and 12 such PDE's for parity-invariant space-periodic flows. 

Let be a solution to the auxiliary problem 



The solvability condition for this equation is orthogonality of the right-hand side to 
the kernel of C For flows symmetric about the vertical axis in a layer, / = 3 and 
satisfies (5) and (1). Solvability is implied by the fact that v is symmetric about 
the axis, and ker£ is spanned by S^ fe ^ + e^, which are antisymmetric about the axis. 
For parity-invariant space-periodic flows, W' ! ' is space-periodic; solvability of (B.l) 
follows from parity-invariance of v and parity-antiinvariance of + spanning 



Thus, it is enough to solve one auxiliary problem (B.l) instead of four problems (14) 
for flows symmetric about the vertical axis in a layer, and three auxiliary problems 
(B.l) instead of nine problems (14) for parity-invariant space-periodic flows. Numer- 
ical complexity of problems (B.l) and (14) is the same, since the operators C and C* 
have the same spectrum. Thus, amount of computations can be reduced twice by the 
use of the auxiliary vector fields W®. 

In computations for parity-invariant space-periodic flows one can further halve 
the number of problems to be solved. The curl of (B.l) is 



£* W (0 



= vxe,. 



(B.l) 



ker C. 



Evidently, from (B.l) and (14) 
(v x T {m ' k) ) l = -(l 




7]V 2 R {1) - V x (v x R«) - = 0, 



(B.2) 



where R w = V x Wj. From (B.2) and (13), 

?|V 2 AW + V x (v x B w ) = 0, 



(B.3) 
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77V 2 B« + V x (v x A«) = -— , (B.4) 

where it is denoted 

A (0 = I(S« + R (/) ), B« = i(S«-R«). (B.5) 

Thus, it is in fact necessary to solve only one equation from the pair (B.3), (B.4), 
e.g. (B.4) in which A^ can be regarded as a notation for r]~ l (— V 2 )~ l V x (v x B^), 
since in the space of Fourier coefficients computation of the inverse Laplacian is a 
simple operation. After and thus A® are determined, subsequently and 
can be found from (B.5), and from (B.l) 

W (0 = (-v 2 )-^ x RW + iy(-V 2 )- 2 V • (v x (R(') + e,)) . 

However, computation of the left-hand side of (B.4) in terms of B^ requires now two 
inverse and two direct FFT's, and thus merging the problems (B.l) and (13) into 
a single equation (B.4) with the notation (B.3) is not necessarily computationally 
advantageous. 

If the flow possesses the "translation antisymmetry" : 

v(x) = — v(x + a) 

for a constant vector a (for the flow (24), (25) this holds for a = ti / L2G2 , if p is odd 
or a 2 = 0), solutions to the problem (13) can be obtained from solutions to (B.l): 

S( fc )(x) = V x W<*>(x + a). 

Alternative representations of elements of the MED tensor are derived in what 
follows using (13) and (14). For j 7^ (and 1 < m, k < 2 for the problem in a layer) 
denote 

£(m,fcj) = e ijx m ( S (k) +Gk + ijr (m ' fe) )/j 2 - 

From (13) and (14) 

C(> (rn,k,3) = jjxm l v l S (k) +f}k + ijT (m,k) + 2 __ \ + ^ x ( y x r (m,fc)) 1 _ ( R g) 

Since the average of horizontal components of the left-hand side of (B.6) vanishes, 
this implies 

V /e ijx ™ L (s« + e k + ^^-) + e m x (v x T^)) ) = 



in 



for any j 7^ 0, and thus 

V J jL + e k + + e m x (v x T^)) dx 3 _ m dx 3 = C (m ' k) A r 



(B.7) 
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Here x 3 - m and x% are Cartesian coordinates in directions orthogonal to x m , integration 
is performed over a section x m =constant of the box of periodicity of the flow, A m is 
the area of the section, and C^ m ' k ^ is a constant vector. Integrating (B.7) in x m and 
normalising by the volume of the box of periodicity, one finds 

= r]e k + V(e rn x (v x T {m ' k) )), 

whereby from (B.7) 

1 cell f)T( m > k ) \ \ 

P(e m x(vxr^)) = — V J J It, f S« + J + e m x (v x T^)j dx 3 . m dx 3 . 

(B.8) 

The right-hand side of (B.8) can be computed for any x m . The m-th component of the 
vector in the right-hand side of (B.8) vanishes as a consequence of (16). Elements of 
the MED tensor (v x r( m,/c )) 3 , required in the problem for a layer, can be determined 
from the only remaining horizontal component of (B.8). 

For a parity-invariant flow, periodic in all three dimensions, the average of the 
left-hand side of (B.6) vanishes, implying by analogy with the case of a layer 

1 r r ( ( f)r( m > k ) \ \ 
(e m x (v x r( m ' fc ))> = —J JL ( S<*> + J + e m x (v x T^) J dx p dx q . 

(B.9) 

Here x p and x q are Cartesian coordinates in directions orthogonal to x m , integration 
is performed over the section x m =constant of the box of periodicity of the flow. The 
right-hand side of (B.9) can be computed for any x m . The m-th component of the 
vector in the right-hand side of (B.9) vanishes as a consequence of (16). The p-ih 
and q-th. components of (v x r*™'^), required in the three-dimensional space-periodic 
problem, can be easily determined from this equation. 

Another expression for (v x r^ m ' fe ^), including the m-th component, can be ob- 
tained from (14) noticing, that vector fields e pq ix q e p (no summation over repeating 
indices is assumed!) solve (B.l) (this observation is, of course, of no consequence for 
determination of W^, since this vector field does not satisfy the required periodicity 
condition (5)). Here e pq i is the permutation symbol: 

G p X G q 6 pq [Qi } 

the three indices p,q and / being distinct. Multiplying (14) by x q , integrating over 
the box of periodicity and normalising by the volume of the box, find 

1 r r ( ( f)r( m > k ) \ \ 

(e q x (v x T^)} = — j j f 77 f 2c^S^ + \ + e,x(vx T^)\ dx t dx p 

+ (x q e m x (vx (S( fc ) + e fc ))>. (B.10) 

Here 5^ is the Kronecker symbol and the integral over a section of the box of pe- 
riodicity of the flow in the left-hand side can be evaluated for any x m . For q = m 
(B.10) does not immediately reduce to (B.9). To perform the reduction, note that 
in the three-dimensional space-periodic problem = 0, cross-premultiply (A. 6) by 
x m e m , integrate over the box of periodicity, normalise the result by the volume of the 
box and subtract it from (B.10). 

Evidently, (B.10) also holds for flows in a layer for 1 < p, q < 2 and I = 3. 
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